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Data from the KTeV experiment have provided many opportunities for studies of fundamental symmetries in na- 
ture, not just CP. We present new Umits on the lepton flavor violating decays — > vr'^/^^eT, — * 



and tt" — >■ ^^e^. We also present measurements of K]^ — > n'^'yy and K]^ — > Tr^e+e 7 which are relevant to 
the determination of the various contributions to Kj^ 



TT^ll. Finally, we present an analysis of the decay 
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1. Introduction 

Symmetries have yielded a significant number of 
insiglits into tire interactions of fundamental parti- 
cles. While CP violation has been intensely studied 
for signs of new physics, only recently has lepton fla- 
vor violation (LFV) exposed cracks in the Standard 
Model. It is not presently clear what kind of mech- 
anism is responsible for the violation of lepton flavor 
symmetry, and the situation is not helped by the fact 
that the violation of this symmetry has thus far ap- 
peared in only the neutrino sector, where experiments 
are extremely challenging. New limits on other mani- 
festations of LFV may be very useful in constraining 
models, while the kaon sector provides a landscape 
where backgrounds to such searches are quite man- 
ageable. However, heavily suppressed CP violating 
decays, such as — ?■ ir'^ll, remain sensitive probes 
of physics beyond the Standard Model owing to the 
ease with which the SM predictions can be produced. 
Finally, CPT remains one of the few conserved sym- 
metries, and as such is an excellent place in which to 
look for new physics. 



and had the effect of transforming the relevant beam 
into a coherent K^-Kg beam. 
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Figure 1: A cartoon of the KTeV detector in the E799 
configuration. The relative positions of the major detector 
components are shown, including the drift chambers and 
Csl calorimeter. Notice the twin kaon beams entering the 
detector. 



3. Search for Lepton Flavor Violation 



,^ 2. The KTeV Experiment 
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KTeV is the culmination of a series of neutral kaon 
experiments at Fermilab. The project consisted of 
two different experiments sharing a common detec- 
tor. E799-II focused on rare Kl decays. E832 was 
designed to establish the existence of Direct CP vi- 
olation in if — >■ TTTT. The KTeV detector, described 
more fully in [H consisted of a long vacuum decay 
volume followed by a magnetic spectrometer, a very 
high performance Csl crystal calorimeter and finally 
a muon veto. KTeV, when operating in E799 config- 
uration, also included a transition radiation detection 
(TRD) system, designed to improve the detector's al- 
ready excellent particle ID capability. In this mode, 
two nearly parallel beams were produced. One the 
other hand, the E832 configuration inserted an active 
regenerator, consisting of a series of plastic scintillator 
blocks instrumented with PMTs, into the beam. This 
regenerator alternated between the two kaon beams. 



With the observation of neutrino oscillation [2] it is 
now known that lepton flavor is not a perfectly con- 
served symmetry of nature. With the root mechanism 
of neutrino oscillation still unknown, it is important 
to search for other manifestations of this symmetry vi- 
olation. Neutral kaon decays can offer an alternative 
probe of LFV. The experimental signature is simple: 
the presence of the lepton pair fi^e^ in any kaon de- 
cay. It should be emphasized that while these modes 
could occur due to neutrino mixing, said rates would 
be undetectable in current experiments Q. Thus ob- 
servation of any new LFV process would provide a 
unique signature of new physics beyond the standard 
model, distinct from neutrino mixing. 



3.1 . Search for Kl n^/j- 



One 



decay mode which would exhibit LFV is 
=F. The E799 data set is filtered [1] to 
look for events in which there is a charged track which 
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matches hits in the muon detector, implying that a 
muon is indeed present. The other charged track 
must then have a Csl cluster whose energy is equal 
to the momentum as measured using magnetic deflec- 
tion, implying that the shower in the Csl was electro- 
magnetic in origin and thus due to an electron. Addi- 
tionally, the transistion radiation detector information 
for the electron-like track must also be consistent with 
an electron. Only events with two charged tracks are 
accepted. Additionally, there must be two more Csl 
clusters not associated with any tracks, which come 
from the two 7s from the tt'^. With these require- 
ments, the signal is dominated by K^i and/or Kf,A de- 
cays where the pion either decays (and thus produces 
an muon) or punches through the Csl calorimeter to 
the muon detector. Making a tight cut on acciden- 
tal activity in the detector eliminates extra Csl clus- 
ters which are needed for Ke'i events to be accepted, 
while Kf-A events can be reduced by calculating |p^|, 
the magnitude of the neutrino momentum assuming 
that the decay was actually a K^-i event, and then 
retaining events for which \p^\ < 0, which is the case 
when the K^A hypothesis is incorrect. 

Using a detailed Monte Carlo simulation of the de- 
cay and the detector, we define contours as a function 
of the Tr'^/ie invariant mass and ( the squared com- 
ponent of the summed momentum which is transverse 
to the neutral kaon beam) which would enclose either 
99% or 95% of the total 7r° /i±e=F sample. The sig- 
nal region is enclosed within the 95% contour and is 
blinded while the control region is defined as falling 
between the 95% and 99% contours. Both contours 
are shown in Figure [2] 

Monte Carlo studies of all possible sources of back- 
grounds indicate that the background expected in the 
control region is 4.21 ± 0.53 events while for the con- 
trol region the expected number is 0.66 ± 0.23 events. 
Unblinding the signal region, we observe no events, 
while we observe 5 events in the control region. This 
leads to a limit on the branching ratio of 

Br{KL 7rV^e=F) < 7.56 x 10"" (1) 

at 90% confidence, which is a factor of 83 lower than 
the previous limit Q. 

3.2. Search for Kl -> yr^TrV^e^F 

Very similar to — >■ ir^ ^r^e"^ is — >■ 7r°7r°/i^e^. 
The additional tt'^ requires the presence of two more 
Csl photon clusters for reconstruction, the side effect 
of which is a strong suppression of the accidental back- 
ground which was an issue for K]^ — > Tr^/i^e^. Due to 
the lower accidental background, it is possible to re- 
lax cuts relative to Kl Tr^yit^e^ in order to raise 
the sensitivity — in this case the tight cut on acciden- 
tal activity in the detector is removed, as is the TRD 
cut on the electron track. After applying these cuts. 
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Figure 2: A plot of the invariant yr^/ie mass versus trans- 
verse momentum squared for the Kl — >■ Tr'^/i^e^ sample. 
The larger (smaller) contour defines the region which en- 
closes 99% (95%) of Kl ^> n^^i^e^ Monte Carlo events. 



and others as detailed in [j| , the primary background 
is Kl Tr'^TT^Tr^^i, tt'^u = 7r° — e"'"e~7 events. In 
order for Kl tt'^tt^tt^d events to contribute to the 
background, one electron must fail to deposit a signifi- 
cant amount of energy into the Csl calorimeter, and at 
the same time, there must be an accidental hit in the 
muon detector which matches the electron track's tra- 
jectory. In order to prevent this, we place a loose re- 
quirement on the TRD information for the muon-like 
track, requiring that the track is not electron-like. Af- 
ter this cut is added, there are an expected 0.44 ±0.23 
background events in the signal region. In the data, 
after all cuts, there are no events in the signal region, 
leading to a branching ratio limit of: 

Br{KL T:°n°fi^e^) < 1.7 x 10"^" (2) 

at 90 % confidence. This is the first limit placed on 
this decay. 

3.3. Search for 7r° ^l'^c^ 

The Kl 7r°7r°/i±e=F search can be easily extended 
by noticing that the decay Kl — > tt^tt^tt'', 7r° — > ^='=e^ 
has the same final state. All that is required to con- 
strain the rate at which this decay occurs is the addi- 
tional requirement that the /^e pair reconstruct near 
the TT*^ invariant mass. After enforcing this additional 
cut, we are led to the branching ratio limit of: 

Br(7r" -» ^±e=F) < 3.59 x 10"^° (3) 

at 90 % confidence. This is a factor of 10 (2) times 
lower than the previous best limits on tt*^ — )■ /i"e^ 
(tt^ A*^e~) @,l3l- Note that this analysis is equally 
sensitive to both charge modes. 
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3.4. Other searches 

While the above neutral kaon decay modes are the 
complete set of LFV modes that KTeV has recently 
investigated ( Kl — >■ //^^^e^e^ was the first LFV de- 
cay studied by KTeV Q ) it should be noted that this 
does not represent all the possible kaon decay modes 
which may exhibit lepton flavor violation. One good 
example is the decay ii^e^ . While KTeV's ex- 

cellent Csl calorimeter resulted in best current limits 
on all LFV modes involving neutral pions, a search 
for this decay was undertaken by the B87f collabo- 
ration using a detector optimized for two body de- 
cays Their resulting branching ratio limit was 
Bi{Kl M=^e=F) < 4.7 X lO^^^ 90% confidence. 

The charged kaon sector also offers excellent probes 

for LFV, such as Rk = r[K±^l.i='l'' \ • "^^^^ quantity 
offers both a very precise standard model prediction 
[lol [ll[ and is sensitive to lepton flavor violating con- 
tributions [H, . Current measurements of Rk have 
uncertainties in the range of 3% while preliminary re- 
suhs from KLOE NA48/2 fTo^ and NA62 (aka 
NA48/3) Jfl] may reach 1%, 2% and 0.5% respectiv- 
ley. Both the KLOE and NA62 collaborations plan 
upgrades to their detectors in an effort to push the 
precision down to the SM uncertainity of 0.04%. 

In addition to searches in the kaon sector, there are 
other areas in which lepton flavor violation may be 
observed. One process which is currently under inves- 
tigation is e^7 ( current limit is Br(/z^ — > 6^7) 
< 1.2 X 10~^^ [31 ) which wiU be probed down to 
the 10~^^ level by the MEG experiment Another 
lepton flavor violating process is the related process 
/i* -f iV — 7> + iV in which the transition involves a 
virtual photon coming from a nucleus instead of an 
emitted real photon. This process, which would have 
been studied by the MECO experiment had the RSVP 
project proceeded, may instead be investigated by the 
Mu2e experiment, an evolution of the MECO design 
which would be sited at Fermilab [l^. That experi- 
ment plans to search for fj,^ + N ^ + N down to 
approximately i?^e < 10""'^^, where i?^e is the ratio of 
the rate of /i* + — > + iV to the rate of muon cap- 
ture for the same given nucleus. This will be a signifi- 
cant improvement over the current limit of 4.3 x lO^^'^, 
obtained using a Ti target [l7j . 

3.5. Final Note 

It should be emphasized that all of the lepton fia- 
vor violating searches detailed here are background 
free. Specifically, the estimated background for 
Kl 7r"^±eT is 0.66, for Kl 7r°7rV^eT the es- 
timate is 0.44 and the expected background for 
7r° ii^e'f is 0.03 events. This implies that a KTeV 
style experiment with additional beam intensity can 
further lower these limits, especially ir^ — >■ /i^e^ due 



to the multiplicity of neutral pions produced via 



4. CP conserving contributions to 

While lepton flavor violation is a direct gateway to 
new physics, various extensions to the Standard Model 
result in a multitude of new or enhanced CP violation 
effects. The suite of decays that comprise Kl tt^U 
are especially attractive targets for further study due 
to the heavily suppressed standard model predictions 
and the precision with which the SM predictions can 
be made. 

However, these modes are also a challenge to ob- 
serve, especially Kl ~> ir'^vv. Kl — >■ 7r°e"'"e~on 
the other hand offers a higher branching ratio and 
a fully reconstructible flnal state, at the cost of a 
less precise standard model prediction. In addition, 
Kl — i> 7r''e'''e~receives contributions from CP con- 
serving, direct CP violating and indirect CP violating 
terms. The challenge is to try to determine the sizes of 
the relative contributions in order to focus on the di- 
rect CP violating component, which is most sensitive 
to new physics effects. Observation of Ks — >■ Tr^e+e" 
constrains the magnitude of the indirect CP violating 
term [18] which then leaves the size of the CP con- 
serving component to be determined. 

4.1. Kl 7r°77 

The vector meson exchange amplitude Ay can be 
used to determine the size of the CP conserving com- 
ponent of Kl — >■ 7r°/Z [l^ In order to produce a sample 
of Kl — ^ 7r°77 events, data from the E832 dataset is 
searched for events containing four Csl clusters, and 
no charged tracks. The energy of each of the Csl clus- 
ters is required to be above 2 GeV. The distribution 
of these photon clusters is also required to be consis- 
tent with a kaon decay within the E832 configuration's 
pure Kl beam. The clusters are then paired such 
that one pair yields an invariant mass within 3 MeV 
of the mass of the tt*^ , while the other pair is required 
to not reconstruct into a tt*^. This last requirement 
suppresses the background due to Kl —> 7r°7r°. Af- 
ter these cuts and others are applied [l^l we are left 
with 1982 events with an expected background level 
of approximately 30%. 

The decay Kl — > tt^tt^ was chosen as the normaliza- 
tion mode, as it has the same final state. The resulting 
branching ratio is: 

Br{KL 7r''77) = (1.29 ± 0.06) x 10"^ (4) 

which is consistent with the result from NA48 .26] and 
supersedes the previous KTeV result j2l| . 
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After having obtained a sample of Kl 7r°77 
the resuhing events were used to perform a likeh- 
hood fit using the model developed in (2^] and the 
Dalitz variables ZjyaUtz 



MIJMk' and Yd 



litz 



{Ej3 — Ery4) / Mk where M34 is the invariant mass of 
the two photons which are not daughters of the tt". 
The result is: 



Av = -0.31 ± 0.05.,t„. ± 0.07, 



syst 



(5) 



The distributions of the Dalitz plot variables for both 
data and Monte Carlo, generated using the best fit 
value of Ay , are shown in Figures |3] and SI 
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Figure 3: Distribution of the variable Yoautz used in the 
analysis of Kl — > 7r°77. 



4.2. Kl ^ 7r°e-' 



e 7 



This decay, which can be though of as n^'-f*j, 
also sees contributions from the vector meson ex- 
change process, and as such can also be used to mea- 
sure Ay. Reconstruction begins by looking for two 
charged tracks and three neutral Csl clusters, which 
are not associated with any tracks. Two of the three 
neutral Csl clusters are required to reconstruct with 
an invariant mass near the mass of the 7r°. The de- 
cay vertex obtained from the tt^ reconstruction is then 
used to compute both M^e-y and Mf-e-y-^-y- In this case 
the decay vertex obtained from the extrapolation of 
the two charged tracks is not used as the two electron 
trajectories are nearly parallel and overlapping in this 
decay. Finally an additional constraint is applied that 
requires that no combination of charged tracks and a 
photon may reconstruct near the invariant mass of the 
tt" . After application of these and additional cuts [23| 



450 
400 
350 
300 
250 
200 
150 
100 
50 




• Data 

itV+ BkgMC 

I I BkgMC 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
m„ GeV/c^ 



Figure 4: Distribution of the variable Znautz used in the 
analysis of Kl — >■ 7r'^77. 



we are left with 139 events with an estimated back- 
ground of approximately 14.4 events. 

Isolating the decay K l ir^ tt^ u i which has the 
same final state, and utilizing it as a normalization 
mode, we then obtain the branching ratio: 



Br{KL Tr°e+e--i) = (1.62 ± 0.17) 



X 10 



-8 



(6) 



which supersedes the older KTeV result [2J] which 
was affected by an older value of the — >■ 7r°7r°£) 
branching ratio. This branching ratio measure- 
ment also compares favorably to the value ob- 
tained using Chiral Perturbation Theory, which is 
Br{KL n°e+e-j)chPT = 1-51 x 10"^ This 
value of the branching ratio also informs searches for 
Kl — TT^e'^e^ , as the low value implies that this de- 
cay mode will not be a sizable contribution to the 
background for Kl — > ir'^e^e^ . Furthermore, the dis- 
tribution of M^Ogg peaks well away from the kaon 
mass in this decay. 

The sample of Kl — > 7r'^e+e^7 events can also 
be used to extract Ay much in the same way as 
Kl — >■ 7r''77. Perfor ming a likelihood fit using the 
model from Reference [221 and the kinematic variables 



^Dalitz 



= M, 



' /Mk' 



ee7 



YLalitz 
2 



and Qoaiitz = MI^/Mk^. Note that Zoaiitz 
and YL)aiitz are analogues of the same variables in 
Kl — > 77^77. The resulting fit produces 



Av 



-0.76 ±0.16,tat± 0.07, 



syst 



(7) 



The distributions of the three kinematic variables, for 
data and Monte Carlo generated with the fitted value 
of Ay, are shown in Figure [S] 
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determine how to form the Dalitz pairs. We choose 
the pairing of the particles in a way such that xi < X2 
and the product xi x X2 is minimized, where xi and 
X2 are kinematic variables defined by: 



(8) 



and are thus proportional to the invariant mass of 
each Dalitz pair. Note that there will be mispairings 
introduced using this method, however this mispairing 
effect is minimized in the full treatment as detailed 
in (soj by utilizing a matrix element method. Once 
the e"'"e~ pairs are formed, the angle is defined as 
the angle between the normal to the plane formed by 



and the plane formed by e 



The resulting 



distribution of the angle <f> is shown in Figure [SI 



Figure 5: Distributions of the three Dalitz variables used 
in the analysis of Kl — > 7r''e^e~7. 

5. The parity of the 7r° 

There has been considerable direct and indirect 
[28l . [29} evidence of the pseudo-scalar nature of the 7r° 
for some time. While the fact that the ir'^ is a pseudo- 
scalar is not in doubt, it is interesting to re-examine 
the issue by studying the decay tt*^ — >■ e'^e~e'^e~ . In 
this decay, each electron-position pair act as an ana- 
lyzer that reveals the polarization of the virtual pho- 
ton from the underlying decay i:^ 7*7*, and thus 
expose the parity of the tt". 

Data from the 4-track trigger of KTeV's E799 
phase were filtered in order to expose events con- 
sistent with the decays Kl — > 'k'^tt^tt^ DD^ t^'^dd = 
tt" e~^e~e~^e~ . The first requirement was that two 
pairs of electromagnetic clusters in the Csl calorime- 
ter, unassociated with charged tracks, each recon- 
struct with an invariant mass consistent with a tt'^. 
The second requirement was that there are also four 
charged tracks in the detector whose energy depo- 
sition in the calorimeter are consistent with elec- 
trons/positrons. The momenta of the four tracks must 
yield a combined invariant mass also consistent with 
a TT^. Finally, the invariant mass of the entire system 
was required to be equal to that of the parent K^. 
The summed momentum of all the daughter particles 
was also required to point back to the kaon production 
target. After these cuts and additional cuts detailed 
in [30| there are 30511 candidate events with an es- 
timated background of 0.6%. The main background 
was from Kl Tr'^j^Tr^/^Tr" where tt'^d — tt''^ ^ e+e~7 

Once the 7r° e'^e~e'^e~ sample has been isolated, 
we can begin to analyze the kinematics of the decay. 
The plane of each Dalitz pair reveals the polarization 
of the parent virtual photon. It is thus necessary to 
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Figure 6: The distribution of angle between the planes de- 
fined by the two Dalitz pairs. There is a strong preference 
for the planes to be orthogonal, consistent with the par- 
ent being a pseudo-scalar particle. The non-zero number 
of events in the regions near = 0, tt, 27r are a result of 
mispairing of the Dalitz pairs. 

In addition to inspecting the distribution of the 
angle (f>, a detailed analysis of the form factor of 
the 7r° — i> e+e^e+e^ decay was carried out. For this 
study, a maximum likelihood fit was carried out us- 
ing the DIP model |31|] and various additional terms 
which allow for scalar and CPT couplings (s^l- The 
fit yields a value for the DIP parameter 

a = 1.3± l.O(stat) ±0.9(sysi) (9) 

which compares favorably to the standard slope pa- 
rameter from 7r° — >■ e+e~7. In addition, the results of 
the fit constrain the presence of a scalar component in 
the decay to < 12.1% (3.3%) at 90% confidence when 
CPT is assumed to be violated (conserved). 

In order to measure the branching ratio of this de- 
cay, we must also select an additional sample in or- 
der to obtain a normalization sample. The decay 
Kl — > TT^ dt^^ dt^^ where tt'^ d — tt'^ ^ e+e~7 is 
selected as a normalization mode due to the simi- 
lar final state. After all cuts, 141251 normalization 
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events remain with a background of 0.5%. The main 
background was Kj^ — >■ ■k^tt'^tt^ . Using the num- 
ber of TT^TT^TT^ events to the number of 
Kl — >■ 7r°£)7r°£)7r° events, along with the previously 
measured value of the Dalitz branching ratio ( which 
is a leading contribution to the systematic error) we 
obtain: 



Brin^ e 



e e+e ) 



(3.26 ± 0.18) X 10-^ (10) 



The above result also includes a detailed treatment of 



radiative corrections [32| . It should be noted that the 
KTeV collaboration intends to release a new measure- 



ment of Br{-K^ 



7) in the near future, at which 



point the above branching ratio can be recomputed 
and the systematic error reduced. 



6. Summary 

The KTeV collaboration has recently produced a 
suite of new limits on the lepton flavor violating decays 
Kl n"fi^e^, Kl iT^TT^ii'^e^ and ^i^e^, 
all of which are either first or current best limits. 
There were no expected backgrounds for all three de- 
cays, indicating that these limits, much like many 
other LFV searches, would be relatively straight for- 
ward to improve in future experiments. In addi- 
tion to LFV searches, the collaboration has also pro- 
duced new measurements of the branching ratios and 
form factors for the related decays Kl — > tt'^jj and 
Kl — > 7r''e+e~7, yielding two new estimates of the 
vector meson exchange amplitude Ay ■ Finally, KTeV 
has produced a new, high statistics analysis of the de- 
cay TT*^ — )> e'^e~e^e~ and its decay characteristics. 
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